ENTROPY-DISSIPATIVE DISCRETIZATION 
OF NONLINEAR DIFFUSION EQUATIONS 
AND DISCRETE BECKNER INEQUALITIES 



CLAIRE CHAIN AIS-HILLAIRET, ANSGAR JUNGEL, AND STEFAN SCHUCHNIGG 

Abstract. The time decay of fully discrete finite- volume approximations of porous- 
medium and fast-diffusion equations with Neumann or periodic boundary conditions is 
proved in the entropy sense. The algebraic or exponential decay rates are computed ex- 
plicitly. In particular, the numerical scheme dissipates all zeroth-order entropies which 
are dissipated by the continuous equation. The proofs arc based on novel continuous and 
discrete generalized Beckner inequalities. Furthermore, the exponential decay of some 
first-order entropies is proved in the continuous and discrete case using systematic inte- 
gration by parts. Numerical experiments in one and two space dimensions illustrate the 
theoretical results and indicate that some restrictions on the parameters seem to be only 
technical. 



1. Introduction 

This paper is concerned with the time decay of fully discrete finite-volume solutions to 
the nonlinear diffusion equation 

(1) Ut = A(w^) in fi, t > 0, u(-,0)=tto in Q, 

and with the relation to discrete generalized Beckner inequalities. Here, (3 > and Q C M, d 
(d > 1) is a bounded domain. When /3 > 1, (TTJ is called the porous-medium equation, 
describing the flow of an isentropic gas through a porous medium [53] • Equation (TTJ with 
/3 < 1 is referred to as the fast-diffusion equation, which appears, for instance, in plasma 
physics with (3 = | [5] or in semiconductor theory with < (3 < 1 [26]. We impose 
homogeneous Neumann boundary conditions 

(2) VfV 3 ) • v = on dQ, t > 0, 
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where v denotes the unit normal exterior vector to dQ, or multiperiodic boundary con- 
ditions (i.e. Q equals the torus T d ). Let us denote by m the Lebesgue measure in M. d or 
jgd-i. we assume for simplicity that m(f2) = 1. For existence and uniqueness results for 
the porous-medium equation in the whole space or under suitable boundary conditions, we 
refer to the monograph [35] . 

In the literature, there exist many numerical schemes for nonlinear diffusion equations 
related to ([1]). Numerical techniques include (mixed) finite-element methods [TJ [151 133] . 
finite- volume approximations [191 E2], high-order relaxation ENO-WENO schemes [TTj . 
or particle methods [30]. In these references, also stability and numerical convergence 
properties are proved. 

The preservation of the structure of diffusion equations is a very important property 
of a numerical scheme. For instance, ideas employed for hyperbolic conservation laws 
were extended to degenerate diffusion equations, like the porous-medium equation, which 
may behave like hyperbolic ones in the regions of degeneracy [31]. Positivity-preserving 
schemes for nonlinear fourth-order equations were thoroughly investigated in the context of 
lubrication- type equations [H|37] and quantum diffusion equations [25]. Entropy-consistent 
finite- volume finite-element schemes for the fourth-order thin-film equation were suggested 
by Griin and Rumpf [22] . For quantum diffusion models, an entropy-dissipative relaxation- 
type finite-difference discretization was investigated by Carrillo et al. [PJ. Furthermore, 
entropy-dissipative schemes for electro-reaction-diffusion systems were derived by Glitzky 
and Gartner [2D]. However, it seems that there does not exist any systematic study on 
entropy-dissipative discretizations for ([I]) and the time decay of their discrete solutions. 

Our first aim is to prove that the finite- volume scheme for ([I])-©, defined in (1301) . 
dissipates the discrete versions of the functionals 



In fact, we will prove (algebraic or exponential) convergence rates at which the discrete 
functionals converge to zero as t — > oo. We call E a a zeroth-order entropy and F a a first- 
order entropy. The functional F\ is known as the Fisher information, used in mathematical 
statistics and information theory [16]. Our analysis of the decay rates of the entropies will 
be guided by the entropy-dissipation method. An essential ingredient of this technique is a 
functional inequality relating the entropy to the entropy dissipation [21 [8] . For the diffusion 
equation (CQ), this relation is realized by the Beckner inequality p[|. 

The entropy-dissipation method was applied to (CQ) in the whole space to prove the decay 
of the solutions to the asymptotic self-similar profile in, e.g., [101 [12]- The convergence 
towards the constant steady state on the one- dimensional torus was proved in [7] . However, 
we are not aware of general entropy decay estimates for solutions to ([1]) to the constant 
steady state, even in the continuous case. The reason might be that generalizations to the 
Beckner inequality, needed to relate the entropy dissipation to the entropy, are missing. As 




(3) 
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our second aim, we prove these generalized Beckner inequalities and provide some decay 
estimates for E a and F a along trajectories of flTj). 

This paper splits into two parts. The first part is concerned with the proof of generalized 
Beckner inequalities and the decay rates for the continuous case. The second — and main — 
part is the "translation" of these results to an implicit Euler finite-volume discretization 
of ([I]). In the following, we summarize our main results. 

The first result is the proof of the generalized Beckner inequality 



q 



(5) Jjf\ q dx- ^Jjfl^dxY <C B (p,q)\\Vf\\ 

where / G and < q < 2, pq > 1. In the case q = 2, we require that \ — -j < p < 1. 

The constant Cb(p, q) > only depends on p, q, and the constant of the Poincare-Wirtinger 
inequality (see Lemma [2] for details). The usual Beckner inequality [I] is recovered for 
q = 2; see Remark [3] for a comparison of related Beckner inequalities in the literature. 
The proof is elementary and only employs the Poincare-Wirtinger inequality. By using a 
discrete version of this inequality (see [6]), the proof can be easily "translated" to derive 
the discrete generalized Beckner inequality 

pq 



\f\ q dx 



r \ pq 

J Q \f\ 1/Pdx ) <C b (p,q)\f\l 



q 

2,T> 



where / is a function which is constant on each cell of the finite- volume triangulation T of 
Q and | ■ | i^,t is the discrete H 1 -seminorm; see Section 13.11 and Lemma [12] for details. 

The second result is the time decay of the entropies E a and F a along trajectories of ([T]). 
Differentiating E a [u(t)] with respect to time and employing the Beckner inequality flSJ), we 
show for > 1 that 

^ [«(*)] < CE a [u(t)]^ + ^ a+1 \ t > 0, 

where C > depends on a, (3, and Cb(p, q)- By a nonlinear Gronwall inequality, this 
implies the algebraic decay of u(t) to equilibrium in the entropy sense; see Theorem [5j If 
the solution is positive and < a < 1, the above inequality becomes 

^fHt)] < C(u )E a [u(t)], t>0, 

which results in an exponential decay rate; see Theorem [61 

The first-order entropies F Q [«(£)] decay exponentially fast (for positive solutions) for all 
(ct, /3) lying in the strip —2 < a — 2(3 < 1 (one-dimensional case) or in the region Ma, 
which is illustrated in Figured] below (multi-dimensional case); see Theorems [7] and |HJ The 
proof is based on systematic integration by parts [22] • In order to avoid boundary integrals 
arising from the iterated integrations by parts, these results are valid only if Q = T d . 
Notice that all these results are new. 

The third — and main — result is the "translation" of the continuous decay rates to the 
finite-volume approximation. We obtain the same results for a discrete version of E a in 
Theorems [141 (algebraic decay) and [T5l (exponential decay). The situation is different for 
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the first-order entropies F a . The reason is that it is very difficult to "translate" the iterated 
integrations by parts to iterated summations by parts since there is no discrete nonlinear 
chain rule. For the zeroth-order entropies, this is done by exploiting the convexity of the 
mapping x t— > x a+l . For the first-order entropies, we employ the concavity of the discrete 
version of dF a /dt with respect to the time approximation parameter. We prove in Theorem 
[16] that for 1 < a < 2 and (3 = a/2, the discrete first-order entropy is monotone (multi- 
dimensional case) and decays exponentially fast (one-dimensional case). We stress the fact 
that this is the first result in the literature on the decay of discrete first-order entropies. 

Throughout this paper, we assume that the solutions to (pQ) are smooth and positive 
such that we can perform all the computations and integrations by parts. In particular, 
we avoid any technicalities due to the degeneracy (f3 > 1) or singularity ((3 < 1) in ([T]). 
Most of our results can be generalized to nonnegative weak solutions by using a suitable 
approximation scheme but details are left to the reader. 

The paper is organized as follows. In Section [2j we investigate the continuous case. 
We prove two novel generalized Beckner inequalities in Section 12. 1[ the algebraic and 
exponential decay of E a [u] in Section T2.2[ and the exponential decay of F a [u] in Section 
12.31 The discrete situation is analyzed in Section [3j After introducing the finite-volume 
scheme in Section 13.11 the algebraic and exponential decay rates for the discrete version of 
E a [u] is shown in Section I3~3"| and the exponential decay of the discrete version of F a [u] 
is proved in Section 13.41 In Section H] we illustrate the theoretical results by numerical 
experiments in one and two space dimensions. They indicate that some of the restrictions 
on the parameters (a, (3) seem to be only technical. In the appendix, a discrete nonlinear 
Gronwall lemma and some auxiliary inequalities are proved. 



It is convenient to analyze first the continuous case before extending the ideas to the 
discrete situation. We prove new convex Sobolev inequalities and algebraic and exponential 
decay rates of the solutions to (fT]). 

2.1. Generalized Beckner inequalities. We assume in this subsection that Q C M. d 
(d > 1) is a bounded domain such that the Poincare-Wirtinger inequality 



for all / G if 1 (f2) holds, where / = m(f2) _1 f n fdx and Cp > only depends on d and 
Q. This is the case if, for instance, Q has the cone property [291 Theorem 8.11] or if dQ is 
locally Lipschitz continuous [361 Theorem 1.3.4]. Suppose that m(f2) = 1 (to shorten the 
proof). Before stating our main result, we prove the following lemma. 

Lemma 1 (Generalized Poincare-Wirtinger inequality). Let < q < 2 and f G if 1 (fi). 



2. The continuous case 



(6) 



||/-/|U 2( o)<Cp||V/|U 2( a) 



Then 
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Proof. Let first 1 < q < 2. The Poincare-Wirtinger inequality fl6]) 

( 8 ) ll/lli2(Q) - ||/||ii(Q) = 11/ - /||l2(h) < c p\\ ^f\\h(n) 
together with the Holder inequality leads to 

(9) ll/ll <c|||v/||| 2(n) + ||/||! 3(n) . 

Here we use the assumption m(f2) = 1. Since q/2 < 1, it follows that 

WfWUn) < {<%\\vf\\bm + ll/llL ( n)) 9/2 < CpWVfWbm + ll/IIW 

which equals ([7j). 

Next, let < q < 1. We claim that 

(10) a 9/2 - a^V -9 / 2 < (a - 6) f//2 for all a > b > 0. 
Indeed, setting c = b/a, this inequality is equivalent to 

1 _ c i-9/2 < (! _ c )9/2 f O rall0<c<l. 

The function #(c) = 1 - c 1 ""/ 2 - (1 - c) 9 / 2 for c £ [0, 1] satisfies g{0) = g(l) = and 
g"(c) = (q/2)(l - q/2){c~ l - q / 2 + (1 - c)^ 2 " 2 ) > for c G (0,1), which implies that 
g(c) < 0, proving Applying (pTO]) to a = ||/||| 2(n) and b = \\f\\ 2 L i {n) and using ©, we 
find that 

(11) 11/1112(0) - 11/11^11/11^ < (WfWhn) - \\f\\l(n)Y /2 < ^l|V/IH 2(n) . 
In order to get rid of the L l norm, we employ the interpolation inequality 

(12) II/IUho) = j[ IfWl^dx < H/lli^ll/ll^n), 

where 6 = q/(2 - q) < 1. Since (2 - q)6 = q and (2 - g)(l - 0) = 2(1 - q), flTI]) becomes 

Il/lll2 (n) - ||/|| 9 i9(n) < C£||V/||k (n)) 
which is the desired inequality. □ 

Lemma 2 (Generalized Beckner inequality I). Let d > 1, < q < 2, pq > 1 or q = 2, 

\ — ^ < P < 1 (0 < p < 1 if d < 2), and let f £ iJ 1 (f2). T/jen t/ie generalized Beckner 
inequality 

pq 



(13) l^dx-^fl^dxY <C B (p,q)\\Vf\\l 2 



(SI) 



holds, where 



C B (p,q) = 2{Pq o 1)C9p zfq<2, C B (p,2) = C 2 P ifq = 2, 
2-q 



and Cp > zs i/ie constant of the Poincare-Wirtinger inequality 
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Remark 3. The case q = 2 corresponds to the usual Beckner inequality [I] 

J \f\ 2 dx - QT \f\ 2/r dxj < C B (p, 2)||V/|| 2 2(n)) 

where 1 < r = 2p < 2. It is shown in [TJ] that the constant Cs{p, 2) can be related to the 
lowest positive eigenvalue of a Schrodinger operator if Q is convex. On the one-dimensonal 
torus, the generalized Beckner inequality (jl~3l for p > and < q < 2 has been derived 
in [7|. In the multi-dimensional situation, the special case p = 2/q was proved in [13]. In 
this work, it was also shown that (1131) with q > 2 and p = 2/q cannot be true unless the 
Lebesgue measure dx is replaced by the Dirac measure. In the limit pq — > 1, (TT31) leads to 
a generalized logarithmic Sobolev inequality (see ffl5|) below). If g = 2 in this limit, the 
usual logarithmic Sobolev inequality [21] is obtained. □ 



Proof of Lemma\^ Let first q = 2. Then the Beckner inequality is a consequence of the 
Poincare-Wirtinger inequality and the Jensen inequality: 

ciiiv/n| 2(n) > ii/ - /n| 2(n) = n/n| 2(n) - n/n| 1(n) > jf / 2 rfx - QT i/i 2 ^ r , 

where 1 — | <r<2(0<r<2if<i<2). The lower bound for r ensures that the 
embedding if 1 (f2) L 2 / r (f2) is continuous. The choice p = r/2 £ [| - yields the 
formulation f fl~3|) . 

Next, let < g < 2. The first part of the proof is inspired by the proof of Proposition 
2.2 in [13]. Taking the logarithm of the interpolation inequality 

ll/IU^n) < 11/11 w^ll/ll^fn) , 
where q < r < 2 and 9{r) = q(2 — r)/(r(2 — q)), gives 

F(r) := -log [ \f\r d x- 9 -^-\og [ \f\^dx-^—^\og[ \f\ 2 dx<0. 
r Jn 1 Jn 1 Jn 

The function F : [q, 2} — > M. is nonpositive, differentiable and F{q) = 0. Therefore, 
F\q) < 0, which equals 

--Aogf \f\"dx+-( f \f\*dx) [ \f\nog\f\dx 



(f Jn Q 



+ 8'(q) (±\og j^f\ 2 dx -Uog jjf\«dx\ <0. 
We multiply this inequality by q 2 J Q \f\ q dx to obtain 

(14) /'| /n ogJ^< i x<-^-||/||l„ (n) log ll/ll ^» 
Jn \\J Wlko) 1 ~ % 



\Li(n) \\J ilz,9(n) 
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Then, we employ Lemma [1] and the inequality log (a; + 1) < x for x > to infer that 

\\J \\LZ{p) \ \\J \\Li(n) J 

Combining this inequality and (fl4|) . we conclude the generalized logarithmic Sobolev in- 
equality 

2 — g 

The generalized Beckner inequality ( Fl3|) is derived by extending slightly the proof of [2 
Corollary 1]. Let 



(15) / \f\no gT ^—dx < _^||V/ llLa(n) . 



G(r) =rlog / |/| 9/r cfe, 



r > 1. 



The function G is twice differentiable with 
G'(r) = 



\f\ q/r dx^j Qf \f\^ r dx\ogJjf\^dx -lJjf\^ r \og\f\dx 



G'\r) = ( L( I \f\^dx) ( / \f\*/'dz I \f\^\\og\f\fdx-( I Ul^logl/M 



r3 V</n / \Jn Jn 



Ix 



The Cauchy-Schwarz inequality shows that G"(r) > 0, i.e., G is convex. Consequently, 
r i — y e G ^ is also convex and r i— >■ if(r) = — (e G ^ — e G ^)/(r — 1) is nonincreasing on 
(1, oo), which implies that 

H{r)<\\mH{t) = -e G ^G\l) = / \f\Hog-j^dx. 

Jn \\j\\L9(p) 

This inequality is equivalent to 

(16) (I ^ - (£ 1 * L ,og 

Combining this inequality and the generalized logarithmic Sobolev inequality (115j) . it fol- 
lows that 

for all < q < 2 and r > 1. Setting p := r/q, this proves ( 1T31) for all pq = r > 1. □ 

For the proof of exponential decay rates, we need the following variant of the Beckner 
inequality. 

Lemma 4 (Generalized Beckner inequality II). Let 0<q<2,pq>l and f e H 1 (Q). 
Then 

(17) ||/||£f n) f / \f\ q dx- ( I \f\ l/p dxY) < C' B (p,q)\\Vf\\l 2{a) , 



where 
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rtPV ~ 1 )°p in < Q<2 
C B M = { 2-q Vl<?<2, 

(pq-l)C 2 P ifO<q<l. 



Proof. By (HM, it holds that for all < q < 2, 



/ l/l 9 logJ^^<^-||/IH 5( n)log 
Jn \\J\\l9(q) 2 ~ <? 



IZI!_^^_9 llfiw i~~ l|/|l ' 2 (") 



I/?(fi) ^ y II J HL9(f2) 

Then, for g > 1, the Poincare-Wirtinger inequality in the version (Q and the inequality 
log (a; + 1) < x for x > yield 

i9(n)l0g i^ - ll/ll -^ log { c2pl jii^ +1 ) - a - ||/|| ^ ||v/|| w 

Taking into account ffTB]) . the conclusion follows for g > 1. 
Let < q < 1. Suppose that the following inequality holds: 

(19) \\f\\l m + ^C|.||V/|| 2 i2(n) - ll/lli^ > 0. 

This implies that, by (fl6|) and for r = pq, 



2 

L 2 (f7) 



2 

1,9(0) 



(pg-i)g„„, fl , ( (2-<i)c 2 P liv/lli > (n 

i«(n) 1Q S „ I, f 1, 2 



<(p?-i)c?.||v/||!. J{n) ||/||«; ( %, 

which shows the desired Beckner inequality. 

It remains to prove ffl9l) . For this, we employ the Poincare-Wirtinger inequality (jSJ) 

C^||V/||i 2(n) > ll/H^n) - ||/||ii (n) 
and the interpolation inequality ()12p in the form 



to obtain 



i«(n) > ll/llK(n)ll/llz?(n) )/fl » e - - 1 » 



L(n) H ~Cp|lWllL2(n) ~ ll/IU 2 (n) 
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We interpret the right-hand side as a function G of Then, setting A = WfW^^) 

G(y)=yV°A^» + 2 -^A- 2 -^>y, 

q q 

G\y) = \ y V^A^'° 2 - q 



0* q 



6 \9 

Therefore, G is a convex function which satisfies G(A) = and G'(A) = 0. This implies 
that G is a nonnegative function on M + and in particular, G(||/||^i, n s) > 0. This proves 
(fl9l) . completing the proof. □ 



2.2. Zeroth-order entropies. Let u be a smooth solution to (JTJ-flSJ) and let u G L°°(Q,), 
inf^tio > in Q. By the maximum principle, < inf^M < u(t) < sup n Mo in Q for t > 0. 
First, we prove algebraic decay rates for E a [u], defined in (j3J). 

Theorem 5 (Polynomial decay for E a ). Let a > and /3 > 1. Let u be a smooth solution 
to CO-© and m e m& inf n w > 0. Then 



where 



(a + l)(a + 0)2 VC7 B (p,g) 

and Cb(p,q) > is the constant in the Beckner inequality for p = (a + /3)/2 and g = 
2(a + l)/(a + /3). 

Proof. We apply Lemma|2]with p = (a + /3)/2 and g = 2(a + l)/(a + /3). The assumptions 
on a and /5 guarantee that < q < 2 and pq > 1. Then, with / = u^ a+l3 ^ 2 , 

«.H ^ (,/ n - +1 - " (W) * ^ (/„ 1^1'^ , ° +1,/, ° +S) 

Now, computing the derivative, 

dt 7n (a + /3) 2 y n ' 1 

< ( 7?^T ) ^] (a+W(a+1) - 



(« + /3) 2 VCb(p,9), 

An integration of this inequality gives the assertion. □ 
Next, we show exponential decay rates. 
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Theorem 6 (Exponential decay for E a ). Let u be a smooth solution to ([I])-®, < a < 1, 
0>O,u o e with inf n u > 0. Then 

E a [u{t)} < E a [u ]e- At , t > 0. 

The constant A is given by 

4a/3 . r , p-i^ 



for (3 > and 



A = — r- f-— -inf(< _1 ) > 0, 



4a/3(g + 1) ^ 
^( P ,g)(a + /3)2 ||Uo|l ^(n)' 



/or/3 > 1. i/ere ; C£(|(a + 1),2) and C' B (p,q) are the constants in the Beckner inequalities 
f fT5j) and (fT7|) . respectively, with p = (a + /?)/2 and g = 2(a + l)/(a + 



Proof. Let /3 > 0. We compute 



-i 



W a+1)/2 | 2 cfe 



dt (a + l) 2 

<--i^-infK-) / |W^V 2 | 2 ^. 
(a + l) 2 n 

By the Beckner inequality ( !T3|) with p = (a + l)/2, q = 2, and / = ■u (cH ~ 1 ) //2 , we find that 

dE a 4a/3 . ^ 

-5r-- CB( P .2)(a + l) '5 f '"° ) £ « 

and Gronwall's lemma proves the claim. The restriction p < 1 in Lemma [2] is equivalent 
to a < 1. 

Next, let /3 > 1. By Lemma H withp = (a+/3)/2, g = 2(a+l)/(a+/3), and / = u( Q+ «/ 2 , 
it follows that 



a+V 

|w||^a+i(Q) I / u a+1 dx — ( / ucfe 



Jq 



Hence, we can estimate 

1/3-1 



dE, 



a 



■ + P) Jn {a + /3y C' B (p,q) 



dt (a 

and Gronwall's lemma gives the conclusion. Note that in the last step of the inequality we 
used that ||u|| LQ +i (n) > ||«||z,i(n) = IKIIz,i(fi)- D 
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2.3. First-order entropies. In this section, we consider the diffusion equation ([1]) on the 
torus Q = T d . We prove the exponential decay for the first-order entropies (j3J. 

Theorem 7 (Exponential decay of F a in one space dimension). Let u be a smooth solution 
to ([!]) on the one-dimensional torus Q = T. Let u G L°°(Q) with inf^no > and let 
a, P > satisfy — 2 < a — 2(3 < 1. Then 

F a [u(t)}<F a [u }e~ At , 0<t<T, 

where 

A = ^ mf («^" 7_1 ) ^f («r a ) > 0, 7 = |(a + - 1), 
where Cp > is £/ie Poincare constant in ((6]). 

Proof. We extend slightly the entropy construction method of [23]. The time derivative of 
the entropy reads as 

f \u«l%{u^u t ) x dx = ~ [ {u^) xx u^-\u% x dx 



dt 2 J n y 



4 

where we introduced 



/ ((I - l) (fl - 1)£ + (f + /? - 2) £ G £, + el) dx, 



SG = ; ?L = ■ 

U U 



This integral is compared to 



((|-i) 2 e G +(7-2)e G ei+ei) 



where, compared to the method of [23], 7 ^ gives an additional degree of freedom in the 
calculations. In the one- dimensional situation, there is only one relevant integration-by- 
parts rule: 

0= !{u a+ ^ul) x dx= [ u a+f, - 1 ((a + P-4:)£ + 3&£ L )dx. 
Jq Jq 

We introduce the polynomials 

(20) s (o = (f-i)o>- m + (f + - 2) + & 

(21) A,(0 = (I - I)' ^ + (7 - 2)&£ L + & 

T(0 = (« + /3-4)^ + 3^l, 

where C, = (£<?, We wish to show that there exist numbers c, 7 6 I (7 ^ 0) and n > 
such that 

5(C) = So(0 + cT(0 - kD (0 > for all f G M 2 . 
The polynomial 5 can be written as S(£) = a\^Q + g^g^l + (1 ~~ where 

ai = -7(7 - 2 ) 2 ^ + (a + /3 - 4)c + -(a - 2)(/3 - 1), 
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a 2 = -(7 - 2)« + 3c + -(a + 2/3 - 4). 

Therefore, the maximal value for k is At = 1. Let k = 1. Then we need to eliminate the 
mixed term £q£l- The solution of a 2 = is given by c = — |(a + 2/3 — 27), which leads to 

If 2, „ A 2 1 



ai = [l ~ g(« + " 1)J - jg(« " 2/3 - l)(a - 2/3 + 2). 

Choosing 7 = |(a + /3 — 1) to maximize a\, we find that «i > and hence S(£) > if and 
only if -2 < a- 2/3 < 1. 

Using the Poincare inequality ([6]) and the maximum principle, we obtain 

= Jvr*e- 1 s {0d X = -?£ J^-\s (o + cT(O)dx 

7^ Qx(0,oo) y Q 

^-^infK+^-^inf^)^. 
Cp n u q 

For the last inequality, we use that (m 7//2 ) x = ^u^~ a ^ 2 (u a ^ 2 ) x , which cancels out the ratio 
a 2 /'-/ 2 . An application of the Gronwall's lemma finishes the proof. □ 

We turn to the multi-dimensional case. 

Theorem 8 (Exponential decay of F a in several space dimensions). Let u be a smooth 
solution to ([Q) on the torus Q = T d . Let w £ L°°(Q) with inf^Mo > and let 

(a, 0) eM d = {(a,/3) e R 2 : (2 - 2a + 2/3 - d + ad) (4 - 4/3 - 2d + ad + 2/3 + 2/3d) > 

and (a - 2/3 - l)(a - 2/3 + 2) < 0} 

(see Figure^). Then there exists A > 0, depending on a, (3, d, u Q , and Q such that 

F a [u(t)}<^[u Q }e~ At , t>0. 

Proof. The time derivative of the first-order entropy becomes 

(22) ^ = -~ / u"/ 2 " 1 A(u a/2 )A(u fS )dx = / vT+^Sodx, 

dt 2 J Q 4 J Q 

where So is defined in (120]) with the (scalar) variables £g — |Vu|/u, = Au/u. We 
compare this integral to 

/ M «+/3-7-l (A(M 7/2 ))2rfx = 2! f u ^-l Dodx , 

Jq 4 
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2 4 6 

Figure 1. Illustration of the set Mj, defined in Theorem [81 for d = 9. 

where Dq is as in ( 1211) and 7 ^ 0. In contrast to the one- dimensional case, we employ two 
integration-by-parts rules: 



= / div (u a+ P- A \Vu\ 2 Vu)dx = [ u a+p - l T x dx, 
Jn Jn 



div(u a+ P- 3 (V 2 u- AI) • Vu)dx = [ 

Jn 



a 



Todx, 



where 



T 1 = (a + (3-4)& + 2Z GHG + e G tL, 

T 2 = (a + (3- 3)£ GHG -(a + /3- 3)$& + & ~ & 

and ^ G hg = h" 3 V« t V 2 hVm, (,h = u _1 ||V 2 w||. Here, ||V 2 w|| denotes the Frobenius norm 
of the hessian. 

In order to compare V 2 m and Am, we employ Lemma 2.1 of [23]: 



, , 1.. , d fVu T V 2 uVu An 

|V^f>-(A^ + — 



Vd 



Therefore, there exists £r 6 R such that 



el 



d 



d d 



GHG 



e 



G 



d 



e L 



R 



+ 



d 
d 



d d 



■el 



where we introduced £s = Cghg/Cg ~%L/d. Rewriting the polynomials T\ and T 2 in terms 
of £=(£g,£l,£s, 6?)eR 4 leads to: 



Ti(0 = (« + /?- 4)£ G + + 2^5, 



1-d 



d 

a + - 3)£ki 



d 



d 



& + ^g(« + - 3) 



R- 



2 

K, 
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We wish to find c 1; C2, 7 6 M (7 / 0) and k > such that 

S(0 = 5 (0 + ciT^O + c 2 T 2 (0 - «D (0 > for all £ e 
The polynomial S 1 can be written as 

S(0 = aid; + °2^g^ + + a ^Us + as£l + c & wnere 

fll = (f _ (/? ~ 1} + {a + P ~ 4)Cl _ (I _ 

01 ( 2 \ , „ , d — 1 

a 2 = - + /3-2+f- + lJ Cl -(a + /3 - 3)^-c 2 - (7 - 2)k, 

1 -d 

a 3 = H — c 2 - ft, 

a 

a 4 = 2ci + (a + /3 - 3)c 2 , 
d 

° 5 = d^l C2 " 

We remove the variable £r by requiring that c 2 > 0. The remaining polyomial can be 
reduced to a quadratic polynomial by setting x = £l/£q arid y = £,s/£g : 

(23) S(x, y) > a,\ + a 2 x + a 3 x 2 + a A y + a 5 y 2 > for all x, y EM.. 

This quadratic decision problem can be solved by employing the computer algebra system 
Mathematica. The result of the command 

Resolve [ForAll [{x, y}, Exists [{CI, C2, kappa, gamma}, 
al + a2*x + a3*x~2 + a4*y + a5*y~2 >= && kappa > 
&& gamma != 0]], Reals] 

gives all (a, (3) E M 2 such that there exist c±, c 2 , 7 E M (7 7^ 0) and n > such that ( 123]) 
holds. This region equals the set M^, defined in the statement of the theorem. 
Similar to the one-dimensional case, we infer that 

dF^^_O^K I u a + p-l Doi0dx = _^ f ^-7-1(^7/2)2^ 

Thus, proceeding as in the proof of Theorem [7] and using the identity 



(Af) 2 dx = / \\V 2 f\\ 2 dx 
Jn 

for smooth functions / (which can be obtained by integration by parts twice), we obtain 



dFa _^/^: ft a+/3- 7 -l\- 

dt S C 2 n [U ° jl n HM ° 



Gronwall's lemma concludes the proof. □ 

Remark 9. Under the additional constraints a 2 = 03 = 0, we are able to solve the decision 
problem (|23|) without the help of the computer algebra system. The solution set, however, 



ENTROPY DISSIPATIVE FINITE- VOLUME SCHEME 15 

is slightly smaller than which is obtained from Mathematica without these constraints. 
Indeed, we can compute c\ and c 2 from the equations a 2 = a 3 = giving 

d fa _ _.\ d(l — k) 

ci = - — - --1 + k(1 + 7-«-/3) , c 2 - 



d + 2 V2 y 1 ' ) d-1 

The decision problem ( 1231) reduces to 

ai + a 4 y + a 5 y 2 > for all y EM.. 

If At < 1, it holds c 2 > and consequently, > 0. Therefore, the above polynomial is 
nonnegative for all y E M if it has no real roots, i.e., if 

< 4aia 5 - a 2 = q + + q 2 'J 2 

for some 7 ^ 0, where (for d > 1) 

d 2 re 



92 



— (3d(rf-4)K+ (d + 2) 2 ) < 0, 



(d + 2) 2 (d 

and go, Qi are polynomials depending on d, a, (3, and k. The above problem is solvable if 
and only if there exist real roots, i.e. if 

0<ql- Aq q 2 = JJ^^y^o + »i« + S 2 K% 

where 

s = -d(5d - 8) + 6d(d - l)a + 2d(d + 2)/3 + 2(d + 2)a/3 - (2d 2 + l)a 2 - (d + 2) 2 /3 2 , 
si = 2d(3d - 4) - 2d(4d - 3)a - 4d(d + l)/3 + 2d(3d - 5) a/3 + 2d(d + l)a 2 

-2d(d-6)/3 2 , 
s 2 = — d 2 (a + /3 — l) 2 . 

We set /(re) = S0 + S1K + S2K 2 . We have to find < re < 1 such that /(re) > 0. Since s 2 < 0, 
this is possible if /(re) possesses two (not necessarily distinct) real roots re and K\ and if 
at least one of these roots is between zero and one. Since /(l) = — (d — l) 2 (a — 2/3) 2 < 0, 
there are only two possibilities for re and K\. either re < < n\ < 1 or < re < k± < 1. 
The first case holds if f(0) = so > 0, the second one if 

(24) f(0) = Sl > 0, /'(l) = Sl + 2s 2 < 0, 

(25) s 2 - 4s s 2 = -4d 2 (a - 2/3 + 2)(a - 2/3 - 1)(4 - 2d + da + 2d/3) 

x (2-d + (d-2)a + 2/3) > 0. 

The set of all (a, /3) G M 2 fulfilling these conditions is illustrated in Figure EJ □ 



3. The discrete case 



We introduce the finite-volume scheme and prove discrete versions of the generalized 
Beckner inequality as well as the discrete decay rates. 
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4 r 



3 - 



2 - 



1 - 




bj i i i i i i i i i i i i i i i lz 

2 4 6 8 

Figure 2. Set of all {a, 0) fulfilling s > 0, and IBS]) for d = 9. 

3.1. Notations and finite- volume scheme. Let Q be an open bounded polyhedral sub- 
set of M. d (d > 2) with Lipschitz boundary and m(fi) = 1. An admissible mesh of Q is 
given by a family T of control volumes (open and convex polyhedral subsets of Q with 
positive measure); a family £ of relatively open parts of hyperplanes in M. d which represent 
the faces of the control volumes; and a family of points (xk)kgt which satisfy Definition 
9.1 in [IT]. This definition implies that the straight line between two neighboring centers 
of cells (xk, %l) is orthogonal to the edge a = K\L between the two control volume K and 
L. For instance, triangular meshes in M. 2 satisfy the admissibility condition if all angles 
of the triangles are smaller than ir/2 [T7J Examples 9.1]. Voronoi meshes in M. d are also 
admissible meshes [TTf Examples 9.2]. 

We distinguish the interior faces a G £mt and the boundary faces a G S ext . Then the 
union £ int U S cxt equals the set of all faces S . For a control volume K G T, we denote by 
£k the set of its faces, by £ lTL t,K the set of its interior faces, and by S ext ,K the set of edges 
of K included in dVt. 

Furthermore, we denote by d the distance in M. d . We assume that the family of meshes 
satisfies the following regularity requirement: There exists £ > such that for all K G T 
and all a G S^k with a = K\L, it holds 

(26) d(x K ,a) > £d(x K ,x L ). 

This hypothesis is needed to apply a discrete Poincare inequality; see Lemma [TH Intro- 
ducing for a G £ the notation 

d _ / d(x K ,x L ) if cr G £ int , o = K\L, 
\ d(x K ,cr) if a G £ cxt ,x, 

we define the transmissibility coefficient 
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The size of the mesh is denoted by Ax = max^gr diam(K) . Let T > be some final time 
and Mt the number of time steps. Then the time step size and the time points are given 
by, respectively, At = j^-, t k = kAt for < k < Mt- We denote by V an admissible 
space-time discretization of Qt = f2 x (0, T) composed of an admissible mesh T of Q and 
the values At and Mt- 

A classical finite volume scheme provides an approximate solution which is constant on 
each cell of the mesh and on each time interval. Let X(T) be the linear space of functions 
Q — y M which are constant on each cell K G T . We define on X(T) the discrete L p norm, 
discrete W 1,p seminorm, and discrete W 1,p norm by, respectively, 

\\u\\o, P ,T= [J \u\ p dx\ = ( ^ m(K)\u K \ 

i i ( m H A 1/P 

<j=K\L 

IM|l,p,T = IM|o,p,T+ |w|l,p,T) 

where u G X(T), u = uk in K G T, and 1 < p < oo. The discrete entropies for u G X(T) 
are defined analogously to the continuous case: 



(27) 



a 



XA'eT VifeT / 



a+l 



(28) F d a [u] = l\u a/2 \l^ r . 

We are now in the position to define the finite- volume scheme of ([I])-®. Let D be a 
finite-volume discretization of Qt- The initial datum is approximated by its L 2 projection 
on control volumes: 

(29) u° = v° k 1k, where u° K = — — — / uo(x)dx, 

and Ik is the characteristic function on K. Then ||w°||o,i,r = H^olU^n)- The numerical 
scheme reads as follows: 

fc+l k 

(30) m(K)^—^ + T.((utty - = 0, 



a=K\L 



for all K G T and /c = 0, . . . , My — 1. This scheme is based on a fully implicit Euler 
discretization in time and a finite- volume approach for the volume variable. The Neumann 
boundary conditions (T5]) are taken into account as the sum in (!30|) applies only on the 
interior edges. The implicit scheme allows us to establish discrete entropy-dissipation 
estimates which would not be possible with an explicit scheme. 



18 C. CHAIN AIS-HILLAIRET, A. JUNGEL, AND S. SCHUCHNIGG 

We summarize in the next proposition the classical results of existence, uniqueness and 
stability of the solution to the finite volume scheme (|29p - (|30p . 

Proposition 10. Let u G L°°(fl), m > 0, M > such that m < u < M in fi. 

Let T be an admissible mesh of Q. Then the scheme (I29p -( l30p admits a unique solution 
(u k K ) K £T,o<k<M T satisfying 

m<u k K < M, for all K G T, < k < M T , 

m( K ) U K = NIU 1 ^)' f° r dl °< k < M T- 

KeT 

We refer, for instance, to [17] and [18] for the proof of this proposition. 

3.2. Discrete generalized Beckner inequalities. The decay properties rely on discrete 
generalized Beckner inequalities which follow from the discrete Poincare-Wirtinger inequal- 
ity [21 Theorem 5]: 

Lemma 11 (Discrete Poincare-Wirtinger inequality). Let Q C R d be an open bounded 
polyhedral set and let T be an admissible mesh satisfying the regularity constraint ([26]). 
Then there exists a constant C v > only depending on d and Q such that for all f G X(T), 

(31) ||/-/||o A r<^l/|i A r 
where f = j n fdx (recall that m(Q) = 1 ) and £ is defined in ( 1261 ). 

Lemma 12 (Discrete generalized Beckner inequality I). Let 0<q<2,pq>l or q = 2 
and < p < I, and f e X(T). Then 

pq 

s n. „\ l f\i 

,2,T 



jf \f\ q dx - U \f\ 1/p dxj 9 < C b (p, q)\f\l 



holds, where 



2(pq - l)Cl CI 
C b (p,q)= ( 2_ )e ,/ 2 */<? <2 > C h (p,2) = -£ ifq = 2. 



C p is the constant in the discrete Poincare-Wirtinger inequality, and £ is defined in fl26l) . 

Proof. The proof follows the lines of the proof of Lemma [2] noting that in the discrete 
(finite-dimensional) setting, we do not need anymore the lower bound on p. Indeed, if 
q = 2, the conclusion results from the discrete Poincare-Wirtinger inequality (l3"Tj) and the 
Jensen inequality. If q < 2, let / G X(7~). Then we have from (fT6l) and (HM 

\f\*dx-( [ \f\^dxY q <(pq-l) [ l/^logJj^-ds 

0,q,T 

g 



/nn\ ^ 2 (pg - 1) \\j\\o,2,T 

( 32 ) — — \\f\\o, q ,T lo £j7W~ 



Q \\T\\o, q ,T 
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We employ the discrete Poincare-Wirtinger inequality (I3ip . 

\\f\\l,2,T - WfWlxT = Wf ~ f\\h,T < C 2 v C X \J\l %T , 
which implies, as in the proof of Lemma [U (see ([9])), that 

ll/ll^r<^r 9/2 |/ll 2 ,r+ll/llo,,,r- 

After inserting this inequality into fl32|) to replace ||/||o,2,r an d using log(x + 1) < x for 
x > 0, the lemma follows. □ 

The following result is proved exactly as in Lemma H] using the discrete Poincare- 
Wirtinger inequality (|3T|) instead of ([6]). 

Lemma 13 (Discrete generalized Beckner inequality II). Let < q < 2, pq > 1, and 
feX{T). Then 

\o2t(I \f\ q dx-( [ \f\ 1/p dx) P9 )<C' b (p, q )\f\l^ r 



holds, where 



C' b (p,q) 



q(pg - l)C 2 p 

^q)T Z/1 " g<2 ' 
(pq - l)C 2 v 

- «yo<<z<i, 



C p is the constant in the discrete Poincare-Wirtinger inequality, and £ is defined in ( 1261) . 

3.3. Zeroth-order entropies. We prove a result which is the discrete analogue of Theo- 
rem Recall that the discrete entropies E^[u k ] are defined in (I27|) . 

Theorem 14 (Polynomial decay of E d ). Let a > and f3 > 1. Let {u k K )K&T,k>o be a 
solution to the finite-volume scheme (1301) with miK & -ru K > 0. Then 

^["^ (ciit + C2) U/. g - 1 » - ^o, 

where 

C2 = £>°]-(/?-l)/(«+l), 

and Cb{p, q) for p = (a + 0)/2 and q = 2(a + l)/(« + 0) is defined in Lemma\1 

Proof. The idea is to "translate" the proof of Theorem [5] to the discrete case. To this end, 
we use the elementary inequality y a+1 — x a+l < (a + l)y a (y — x), which follows from the 
convexity of the mapping x i— > x a+1 , and scheme fl30l) : 

E d a [u k+1 ] - E d a [u k ] ~E m(JO((4 +1 ) Q+1 - 
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u k K +1 T(u k K +1 -u k K ) 



<E#)(' 

KeT 

<~ At H E ^(4 +1 ) a ((4 +1 )'-K +1 )0- 



KGTaGS int , 
a=K\L 



Rearranging the sum leads to 



(33) E d a [u k+1 ] - E d a [u k ] < -At £ - K +1 r) - 



a=K\L 



Then, employing the inequality in Lemma [19] (see the appendix), it follows that 
K[u k+l ]-Ei[u k ]<-^^ 2 r ff ((^)W)/ 2 -(^)^/ 2 ) 2 



a=K\L 



< |( M fc+l)(a+^)/2|2 



(a + /3) 2 '^ ' W 
and applying Lemma [U with p = (a + P)/2, q = 2(a + \)/{a + /3), and / = («+0)/ 2 , 

The discrete nonlinear Gronwall lemma (see Corollary [IS] in the appendix) with 

_ Aa(5At ( a + 1 \ _ a^f3_ 

T ~ (a + /3) 2 \C b (p,q)J ' 7 " a + 1 > ' 

implies that 

F d uM < - h > n 

o[ J " (£>T~ 7 + c 1 t fe ) 1 /(7-D ' 

where ci = (7 — 1)/(1 + 7T.E^[m ] 7_1 ). Finally, computing ci shows the result. □ 

The discrete analogue to Theorem [6] is as follows. 

Theorem 15 (Exponential decay for E d ). Let {u k K )K<^T,k>o be a solution to the finite- 
volume scheme ([30]) and let < a < 1, f3 > 0, mfKeT u< k — 0- Then 

E d a [u k ] < Ei[u°}e~ xt \ k>0. 

The constant A is given by 

Aa /3 



A 



C 6 (i(a + l),2)(a+l) 



- inf ((O^ 1 ) > 0, 



/or /? > ; and 

_ 4a/3(a + 1) ^ 
= ^(p,g)(a + /3) 2l|n ll0 ^ 
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for (3 > 1. i/ere C' b (p,q) > is the constant from Lemma [731 wift p = (a + (3)/2 and 
q = 2{a + l)/{a + (3). 

Proof. Let a < 1 and /? > 0. As in the proof of Theorem [TU we find that (see f]33|) ) 

- ^[«*] < -At ^((4 +1 r-(4 +1 r)((4 + Y-(^ + Y)- 

CT=FS:|L 

Employing Corollary [20] (see the appendix), we obtain 

Et[u k ^-Ei[u k ]<-^^ 2 £ ^minK^Y-S^Y- 1 } 

o-=ft:|L 

x ((^ +1 )( Q+1 )/ 2 -(4 +1 )( Q+1 )/ 2 ) 2 
4a/3At fc+i^-ii/ fc+iN(a+i)/2i2 

< 4a(3At , fe+1 

C&^a + l),2)(a + 1) KeT 

where we have used Lemma [T2l with p = (a + l)/2, q = 2, and / = n^ a+1 ^ 2 . Now, the 
Gronwall lemma shows the claim. 

Next, let /3 > 1. As in the proof of Theorem [14], we find that 

Ei[u^] - E d a [u k ] < -^^M +1 ) {a+m \h,r- 

We apply Lemma [Lj with p = {a + 0)/2, q = 2(a + I) /{a + and / = u^)/ 2 to obtain 

E a [u }-E a [u}< E a [u } 

4a/3(a+l)Atll« llfc 1 r d k+1 
(a + py C' b (p,q) a[ J - 

Then Gronwall's lemma finishes the proof. □ 

3.4. First-order entropies. We consider the diffusion equation (TjQ) on the half open unit 
cube [0, l) d C M. d with multiperiodic boundary conditions (this is topologically equivalent 
to the torus T d ). By identifying "opposite" faces on dQ, we can construct a family of 
control volumes and a family of edges in such a way that every face is an interior face. 
Then cells with such identified faces are neighboring cells. 

Theorem 16 (Exponential decay of F d ). Let {u k K )K^T,k>o be a solution to the finite-volume 
scheme (|30~|) with f2 = T d and ia£KeT u< k — 0- Then, for all 1 < a < 2 and f3 = a/2, 

Ft[u k+1 ] < F d [u% keN. 
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Furthermore, if d = 1 and the grid is uniform with N subintervals, 

F d a [u k ] < F d a [u ]e- Xt \ 
where X = 4/3 sin 2 (vr/iV) min^^) 2 ^ 1 )) > 0. 
Proof. The difference G a = F^[u k+1 ] — F^[u k ] can be written as 

o a = \ E ((« +1 ) a/2 - K +1 r /2 ) 2 - ((u k K y /2 - (4r /2 ) : 

<T=K\L 

Introducing clk = {u^ 1 — u k K )/T, we find that 

Ga = \ E To (((4 +1 )° /2 " (utr /2 ) 2 ~ ({u k K +1 ~ ra K r /2 - K+ 1 - ra L )^y 

a=K\L 

We claim that G a is concave with respect to r. Indeed, we compute 

o-=i<:|L 

x - ra^^V - K +1 - ra L ) Q / 2 ^a L ), 



<9 2 G Q a 
Th 2 " = ~1 

cr=K|Z, 



Y E ^((«^ 1 -™Jf) a/a "W-(ti^ 1 -ra i )°^- 1 aJ a 



(T=K\L 

x ((n^ 1 - ™ x ) a / 2 - 2 a 2 , - (u k L +1 - ra L T l2 - 2 a\). 

Replacing u k ^ 1 — rax, u*l +1 — tol by u k K) u\, respectively, the second derivative becomes 

d 2 G a a 2 sr^ , k , a / 2 _i /„,fc\a/2-i„ \ 2 
-q^T = ~ 1^ t °\\ u k) a K -{u L ) a L ) 

°"££int; 

a=K\L 

-f (f- 1 ) E ^((^) a/2 -(^) Q/2 )((^) a/2 - 2 4-K) a/2 - 2 a!) 

<r=ir|z 

= _ T E r ^( c i a i<: + c 2«^ a L + c 3 a 2 ), 
o-=js:|i 



where 



ci = (a - 2)((^)«/ 2 - («|) a/2 )(^) a/2_2 + «(4)°- 2 , 

fe \a/2-l/,.fcW2-l 



c 2 = -2a«)^-> 
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c 3 = -(a-2)[(u K ) 1 - K) 1 ){u L ) 1 +a[u L ) . 

We show that the quadratic polynomial in the variables and is nonnegative for all 
u h K and u k L . This is the case if and only if c\ > and 4ciC3 — c 2 > 0. The former condition 
is equivalent to 

2(a - l)(4) a ~ 2 > (a - 2)(4) a / 2 - 2 (4)«/ 2 , 
which is true for 1 < a < 2. After an elementary computation, the latter condition 
becomes 

4 Cl c 3 -c 2 2 = 8(a - 1)(2 - a)«) Q/2 ^ 2 (4) a/2 " 2 ((^) a/2 - (u k L ) a/2 ) 2 > 

for 1 < a < 2. This proves the concavity of G tt (r). 
A Taylor expansion and G a (0) = leads to 

G a (r)<G a (0) + r^(0) 
or 

= T ^ r -((^ 1 ) Q/2 " « +1 ) Q/2 ) (« 1 ) a/2 - 1 a«: " « +1 r /2 "V) 

= jE ^(4 +1 ) Q/2 -K +1 ) a/2 )(4 +1 ) a/2 - 1 «, 

o-=J<r|i 

+ yE r CT ((4 +1 ) a/2 -(4 +1 ) Q/2 )(4 +1 ) a/2 - 1 «x. 

ce£mt > 
o-=ftr|i 

Replacing and by scheme (I30I) and performing a summation by parts, we infer that 

T 



(34) x («y - «+y ) (« +l r /2 - (4 +1 r /2 ) ■ 

Note that the expression on the right-hand side is the discrete counterpart of the integral 



u^-\u^) xx (u a / 2 ) xx dx, 



appearing in (1221) . The condition a = 2(3 implies immediately the monotonicity of k i— > 
F d a [u k }. 

For the proof of the second statement, let d — 1 and decompose the interval Q in N 
subintervals K t , . . . ,K N of length h > 0. Because of the periodic boundary conditions, 
we may set u k N+l = Uq and u_ x = u k N) where u k is the approximation of the mean value 
of u(-,t k ) on the subinterval K i: i = 1,...,N. We rewrite (134]) for a = 2(3 in one space 
dimension: 

G ^)<-^£( £ - ( u j k+1 f)X 

i=1 ie{i— i,i+i> ' 
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i=l 

where Z{ = (u k+1 Y — (u k ^lY . The periodic boundary conditions imply that J2iLi z i = 0- 
Hence, we can employ the discrete Wirtinger inequality in [3U Theorem 1] to obtain 

N 



^sin 2 |- min ((^) 2(/3 " 1) ) E 
h N i=l,...,N v ; ' 



i=l 



4/3r . 2 vr 



h ' N i=i,...,N 



((u*r^)F«[u k+i ] 



By the discrete maximum principle, m.axi(u k+1 ) 2<yl ^ < max^-u?) 2 *- 1 ^ which is equivalent 
to minj(-u( s+1 ) /3_1 > minj(u°)^ _1 - Therefore, 

F d > k+l ] ~ F d a [u k ] = G 2 ,(At) < sin 2 £ ({u° i ) 2 ^)F^[u k+1 }, 

and Gronwall's lemma finishes the proof. □ 

4. Numerical experiments 

We illustrate the time decay of the solutions to the discretized porous-medium ((3 = 2) 
and fast- diffusion equation ((3 = 1/2) in one and two space dimensions. 
First, let (3 = 2. We recall that the Barenblatt profile 

ti (rt)-(t + t)- A (c \x-Xo\ 2 \ VP" 1 ? 

is a special solution to the porous-medium equation in the whole space. (Here, z + denotes 
the positive part of a function z + := max{0, z}.) The constants are given by 

A=* . B 1 



d(/3-l) + 2' d(/3-l) + 2' 

and C is typically determined by the initial datum via J n u(x,t)dx = f n u(x,0)dx. We 
choose C = B((3 — l)(2/3) _1 (t! + t )~ 2B \xi — x \ 2 , where t\ > is the smallest time for 
which u(xi, ti) = 0. 

In the one-dimensional situation, we choose Q = (0, 1) with homogeneous Neumann 
boundary conditions and a uniform grid (xi,P) G [0,1] x [0,0.2] with 1 < i < 50 and 
< j < 1000. The initial datum is given by the Barenblatt profile ub{-, 0) with xq = 0.5, 
x\ = 1 and to = 0.01. The constant C is computed by using t\ = 0.1, which yields 
C ~ 0.091. For < t < 0.1, the analytical solution corresponds to the Barenblatt profile. 

The time decay of the zeroth- and first-order entropies are depicted in Figure [3] in semi- 
logarithmic scale for various values of a. The decay rates are exponential for sufficiently 
large times, even for a > 1 (compare to Theorem [15!) and for a ^ 2(3 (see Theorem [T6l) , 
which indicates that the conditions imposed in these theorems are technical. For small 
times, the decay seems to be faster than the decay in the large-time regime. This fact 
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has been already observed in [TJ Remark 4]. There is a significant change in the decay 
rate of the first-order entropies F£ for times around t± = 0.1. Indeed, the positive part of 
the discrete solution, which approximates the Barenblatt profile ub for t < t±, arrives the 
boundary and does not approximate ub anymore. The change is more apparent for a < 1. 




-20 



0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 

Figure 3. The natural logarithm of the entropies \og(E^[u](t)) (left) and 
\og(F^[u](t)) (right) versus time for different values of a (/9 = 2, d — 1). 

Next, we investigate the two-dimensional situation (still with /3 = 2). The domain 
Q = (0, l) 2 is divided into 144 quadratic cells each of which consists of four control volumes 
(see Figure H]). Again we employ the Barenblatt profile as the initial datum, choosing 
t = 0.01, ti = 0.1, and x = (0.5,0.5), and impose homogeneous boundary conditions. 
The time step size equals At = 8 • 10 -4 . 




Figure 4. Four of the 144 cells used for the two-dimensional finite- volume scheme. 

In Figure El the time evolution of the (logarithmic) zeroth- and first-order entropies are 
presented. Again, the decay seems to be exponential for large times, even for values of a 
not covered by the theoretical results. At time t = t\, the profile reaches the boundary 
of the domain. In contrast to the one-dimensional situation, since the radially symmetric 
profile does not reach the boundary everywhere at the same time, the time decay rate of 
does not change as distinct as in Figure [31 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 

Figure 5. The natural logarithm of the entropies log(E^[u](t)) (left) and 
log(F^[u](t)) (right) versus time for different values of a (/? = 2, d — 2). 

Let f3 = 1/2. The one- dimensional interval Q = (0, 1) is discretized as before using 51 
grid points and the time step size is At = 2 • 10 -4 . We impose homogeneous Neumann 
boundary conditions. As initial datum, we choose the following truncated polynomial 
u (x) = C((x — x)(x — £i))+, where x = 0.3, X\ = 0.7, and C = 3000. In the two- 
dimensional box Q = (0, l) 2 , we employ the discretization described above and the initial 
datum u (x) = C(R 2 — \x — x | 2 )+, where R = 0.2, x = (0.5, 0.5) and again C = 3000. 

In the fast- diffusion case < 1, we do not expect significant changes in the decay rate 
since the initial values propagate with infinite speed. This expectation is supported by 
the numerical results presented in Figures M and [7J For a large range of values of a, the 
decay rate is exponential, at least for large times. Interestingly, the rate seems to approach 
almost the same value for a G {0.5, 1, 2} in Figure t3 



Appendix A. Some technical lemmas 

A.l. Discrete Gronwall lemmas. First, we prove a rather general discrete nonlinear 
Gronwall lemma. 

Lemma 17 (Discrete nonlinear Gronwall lemma). Let f G C 1 ([0,oo)) be a positive, non- 
decreasing, and convex function such that 1/f is locally integrable. Define 

r dz 

Let (x n ) be a sequence of nonnegative numbers such that x n+1 — x n + f(x n+ i) < for 
n G No- Then 

x n < w- 1 (w(x ) - Yq^y) , neK 



ENTROPY DISSIPATIVE FINITE- VOLUME SCHEME 




Figure 6. The natural logarithm of the entropies log(E^[u](t)) (left) and 
\og(F^[u]{t)) (right) versus time for different values of a (/? = 1/2, d — 1). 




0.02 0.04 0.06 0.08 0.1 




0.02 0.04 0.06 



Figure 7. The natural logarithm of the entropies \og(E^[u](t)) (left) and 
\og(F^[u](t)) (right) versus time for different values of a ((3 = 1/2, d = 2). 



Notice that the function w is strictly increasing such that its inverse is well defined. 
Proof. Since / is nondecreasing and (x n ) is nonincreasing, we obtain 

w(x n+ i) - w(x n ) = 



+ dz x n _|_i x n 



f(z) - /(*„) 
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The sequence (x n ) satisfies f(x n+ i)/(x n+ i — x n ) > — 1. Therefore, 

/ \ / s ^ ( f&n+i) , f(x n ) - /(x n+ i)\-i 
w(x n+ i) - w{x n ) < ( 1 J 

fM - /(X n+ i)\-! 



< 

•^n ^n+l 

By the convexity of /, f(x n ) - f(x n+1 ) < f(x n )(x n - x n+1 ) < f(x )(x n - x n+1 ), which 
implies that 

w{x n+1 ) - w(x n ) < (-1 - /'(xo))^ 1 . 
Summing this inequality from n — to N — 1, where iVeN, yields 

Applying the inverse function of w shows the lemma. □ 

The choice f(x) = tKx 1 for some 7 > 1 in Lemma [T71 lead to the following result. 
Corollary 18. Let (x n ) be a sequence of nonnegative numbers satisfying 

x n+ i -x n + rxl +1 < 0, n E N, 
where K > and 7 > 1 . Then 

Xn - 77^ — 1 x 1/(7-1) ' n e N ' 

where c = (7 — 1)/(1 + 7TXo _1 ). 

A. 2. Some inequalities. We show some inequalities in two variables. 
Lemma 19. Let a, (3 > 0. Then, for all x, y > 0, 

(35) (y a - x a ){y? - x?) > ^ {y {a+m - x^' 2 ) 2 . 

( I p J 

Proof. If y = 0, inequality (j35p holds. Let y 7^ and set z = (x/y)@. Then the inequality 
is proved if for all z > 0, 

f(z) = (1 - ^)(1 - 2 ) - ^^(1 - z {a+m "? > 0. 



We differentiate / twice: 

4a 

_ a{a - 0) g/20-3/2 ( 1 .a/2/3-1/2 , a ~ ^ .a/2/3+1/2 , 2 \ 
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Then /(l) = and /'(l) = 0. Thus, if we show that / is convex, the assertion follows. In 
order to prove the convexity of /, we define 

n ( 7 \ - _iw^-l/J i 01 ~ & r a/2/3+l/2 , 2 
9[) ~ (3 + /?(c* + /3) + 

Then g(l) = and it holds 

and therefore, g'(l) = 0. Now, if a > /?, g(0) = 2/ (a + f$) > 0, and g is decreasing in [0, 1] 
and increasing in [1, oo). Thus, g(z) > for all z > 0. If a < (3 then g(0+) = — oo, and g 
is increasing in [0, 1] and decreasing in [1, oo). Hence, g(z) < for z > 0. Independently 
of the sign of a — (3, we obtain 

/"(*) = a{a ~ P) z a l 2 ^l 2 g(z) > 

for all z > 0, which shows the convexity of /. □ 
Corollary 20. Let a, (3 > 0. T/ien, /or all x, y > 0, 

> min{^- 1 ,/- 1 }(y( a+1 )/ 2 - ^ a+1 >/ 2 ) 2 . 

a + lr 



Proof. We assume without restriction that y > x. Then we apply Lemma [T9l to /3 — 1: 

(/ - ^)(^ - x a ) = ^^(y" - - x) > 4a - J (a+m _ x ( Q+ i )/2)2 _ 

?/ — x (a + 1)^ y — X 

/3 j V t^dt > f3Ynm{xP-\yP- l }(y-x), 

J X 

the conclusion follows. □ 



Since 



yP-x 13 
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